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Abstract. The maximum possible number of non-overlapping unit spheres 
that can touch a unit sphere in n dimensions is called kissing number. The 
problem for finding kissing numbers is closely connected to the more general 
problems of finding bounds for spherical codes and sphere packings. We 
survey old and recent results on the kissing numbers keeping the generality 
of spherical codes. 


1 Introduction 

How many equal billiard balls can touch (kiss) simultaneously another billiard 
ball of the same size? This was the subject of a famous dispute between Newton 
and Gregory in 1694. The more general problem in n dimensions, how many non¬ 
overlapping spheres of radius 1 can simultaneously touch the unit sphere 
is called the kissing number problem. The answer is called kissing number, 
also Newton number (in fact, Newton was right, without proof indeed, with his 
answer T 3 = 12 ) or contact number. 

Further generalization of the problem leads to investigation of spherical codes. 
A spherical code is a non-empty finite subset of Important parameters of a 

spherical code C C are its cardinality \C\, the dimension n (it is convenient 
to assume that the vectors of C span R"") and the maximal inner product 

s(C) = max{(x,y) : x, y G G, x 7 ^ y}. 

The function 

A(n, s) = max{|G| : C with s(G) < s} 
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extends r^, and it is easy to see that A{n,1/2) = r„. One also considers the 
function 

D{n,M) = max{d(C) = 2{1 - s{C)) : 3C C with \C\ = M] 

which is used in the information theory (cf. |13l [TSl |32] ). 

For n > 3 and s > 0, only a few values of A{n, s) are known. In particular, 
only six kissing numbers are known: ri = 2, r 2 = 6 (these two are trivial), 
Ts = 12 (some incomplete proofs appeared in 19th century and Schiitte and van 
der Waerden m first gave a detailed proof in 1953, see also [271 IMl El EZ]), 
r 4 = 24 (hnally proved in 2003 by Musin [37]), rg = 240 and r 24 = 196560 (found 
independently in 1979 by Levenshtein [30] and Odlyzko-Sloane |38jl. 

Note that Kabatiansky and Levenshtein have found an asymptotic upper 


bound for Tn [25] • (Currently known the lower bound is 

m) 


This survey deals with the above-mentioned values of Tn and mainly with 
upper and lower bounds in dimensions n < 32. Some interesting advances during 
the last years are described. 

Usually the lower bounds are obtained by constructions and the upper bounds 
are due to the so-called linear programming techniques and their extensions. We 
describe constructions which often lead to the best known lower bounds. The 
upper bounds are based on the so-called linear programming [H [25] and its 
strengthening [371E91137]. Applications were proposed by Odlyzko-Sloane [38] , 
the first author [7], and strengthening by the third author m and Pfender [39]. 

Recently, the linear programming approach was strengthened as the so-called 
semi-definite programming method was proposed by Bachoc-Vallentin [5] with 
further applications by Mittelmann-Vallentin |33j . 

2 Upper bounds on kissing numbers 

2.1 The Fejes Toth bound and Coxeter-Bdrdczky bonnd 

Fejes Toth [22] proved a general upper bound on the minimum distance of a 
spherical code of given dimension and cardinality. In our notations, the Fejes 
Toth bound states that 



( 1 ) 
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where ipM = 6( m -2) • This bound is attained for M = 3,4,6, and 12. This gives 
four exact values of the function D{n,M) (but not necessarily implying exact 
values for A{n,s)). 

First general upper bounds on the kissing numbers were proposed by Coxeter 
|14j and were based on a conjecture that was proved later by Boroczky [ 6 ]. Thus 
it is convenient to call this bound the Coxeter-Boroczky bound. 

Let the function Fn{a) be defined as follows: 

Fo(a) = Fi{a) = 1, 

Fn+i{a) = - [ Fn-i{P{t))dt 

./(1/2) arccos(l/n) 

for n > 1, where /3(t) = ^ arccos i_^ 2 cos 2 r This function was introduced by 
Schlafli [lO] and is usually referred to as Schlafli function. 

In terms of the Schlafli function the Coxeter-Boroczky bound is 

2 F -\ ( 

T(n, s) < AcB{n, s) = , (2) 

Fn{a) 

where a = ^ arccos ■ 

The bounds < Acsin, 1/2) are weaker than the linear programming bound 
to be discussed below. On the other hand, we have 

A{4, cos 7 r/ 5 ) = 120 = Acs (4, cos vr/S) = 

T4(vr/5) 

(the lower bound is ensured by the 600-cell). The value A(4, cosvr/S) = 120 can 
be found by linear programming as well [2]. This suggests that the Coxeter- 
Boroczky bound can be better than the linear programming bounds when s is 
close to 1 . 

2.2 Pure linear programming bounds 

The linear programming method for obtaining bounds for spherical codes was 
built in analogy with its counter-part for codes over finite fields which was de¬ 
veloped by Delsarte m- Delsarte-Goethals-Seidel [TU] proved in 1977 the main 
theorem and it was generalized by Kabatianskii-Levenshtein [25] in 1978. 

The Gegenbauer polynomials mm play important role in the linear program- 

(n) 

ming. For fixed dimension n, they can be defined by the recurrence Pq ' = 1, 
= t and 

{k + n- 2)P^"^^{t) = {2k+ n- 2)tPlf'\t) - A:P^/;\(t) for A: > 1. 
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If 

m 

fit) = aif 

i=0 

is a real polynomial, then /(t) can be uniquely expanded in terms of the Gegen- 
bauer polynomials as 

m 

f{t) = YfkPt\t)- 

k=0 

The coefficients fi, i = 0,1,..., fc, are important in the linear programming 
theorems. 


Theorem 1. (Delsarte-Goethals-Seidel [TB], Kabatianskii-Levenshtein |25| 1 
Let f{t) be a real polynomial such that 
(Al) fit) < 0 for —l<t<s, 

(A2) The coefficients in the Gegenbauer expansion fit) = YfJk=o fkP^\t) 
satisfy fo > 0, fk > 0 for i = 1,... ,m. 

Then A(n, s) < /(l)//o- 




There are two cases, in dimensions eight and twenty-four, where only techni¬ 
calities remain after Theorem 1. The lower bounds rg > 240 and r 24 > 196560 
are obtained by classical configurations and the upper bounds are obtained by 
the polynomials 

h 

and 

= (t + l)(t + l)2(t + 1)V(t 

respectively (the notations will become clear later). Indeed, one may easily 
check that these two polynomial satisfy the conditions (Al) and (A2) for the 

corresponding values of n and s and therefore rg < = 240 and r 24 < 

/r^_ 


fo 


= 196560. 


Together with the Gegenbauer polynomials we consider their adjacent poly¬ 
nomials which are Jacobi polynomials P^’^\t) with parameters 

, , n — 3 , n — 3, 

(a, /3) = (a H--—, b -\ --—) 


where a,b G {0,1} (the Gegenbauer polynomials are obtained for a 
Denote by the greatest zero of the polynomial P'^’^\t). Then 


t 


1,1 

fc-i 


< t 


1,0 


< t 


1,1 

k 


b = 0). 
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for every k >2. 
Denote 


T — 
-^m — 


1,1 1,0 
''fc-1’ ''fc 
, 1,0 , 1,1 


if m 
if m 


2k - 1, 
2k, 


for k = 1 , 2 ,... and Xq = [— 1 , 

Then the intervals Im are consecutive and non-overlapping. For every s G Im, 
the polynomial 


r(n,s) 
J m 


it) 


{t - s) (rl’^it, , if m = 2 A: - 1 , 

+ s)(Tl']:^{t,s)^ , ifm = 2k, 


can be used in Theorem 1 for obtaining a linear programming bound. Levenshtein 
[30] proved that the polynomials satisfy the conditions (Al) and (A2) 

for all s G Im- Moreover, all coefficients fi, 0 < i < m, in the Gegenbauer 
expansion of fm’^\t) are strictly positive for s G Im- Hence this implies (after 
some calculations) the following universal bound. 

Theorem 2. (Levenshtein bound for spherical codes |30l I31j l Let n > 3 and 
s G [—1,1). Then 


A{n, s) < < 


L2k-iin,s) = ( 


fc+n—3 
fc -1 


) 


2fc+n—3 
72—1 




for s G l 2 k-i, 


L2kin,s) = C+l-^) 


2k-\-n—l 

72—1 


{l-s)lpt'>is)+Pil\{s)) 

for s £l 2 k- 


In particular, one has Tg < 1^(8, 1/2) = L 7 ( 8 ,1/2) = 240 and T24 < Lio(24,1/2) 
Lii(24,1/2) = 196560. The Levenshtein bound can be attained in some other 
cases (cf. the tables in [Sni EH E2] ) • 

The possibilities for existence of codes attaining the bounds Lmin,s) were 
discussed in m- In particular, it was proved in ini Theorem 2.2] that the even 
bounds L 2 k{n,s) can be only attained when s = or s = This follows 
from a close investigation of the two-point distance distribution 


A* = 




^ |{y G C : (x,y} = t}| = |^|{(x,y) G : {x,y) = t}| 


xdC 
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of the possible (n, L 2 k{n, s), s) codes. 

On the other hand it was proved by Sidelnikov [52] (see also |32l Theorem 
5.39]) that the Levenshtein bounds are the best possible pure linear programming 
bound provided the degree of the improving polynomial is at most m. This 
restriction was later extended by Boyvalenkov-Danev-Boumova [9| to m + 2 and 
the polynomials are still the best. 

However, in some cases the Leveshtein bounds are not the best possible pure 
linear programming bounds. This was firstly demonstrated in 1979 for the kissing 
numbers by Odlyzko-Sloane [38]. Boyvalenkov-Danev-Boumova [9] proved in 1996 
necessary and sufficient conditions for existence of certain improvements. 

Theorem 3. [9| The bound Lm,{n,s) can he improved by a polynomial from 
^n,s of degree at least m + 1 if and only if Qj(n, s) < 0 for some j > m + 1. 
Moreover, if Qj{n, s) < 0 for some j > m + 1, then Lm{n, s) ean be improved by 
a polynomial from An,s of degree j. 

For s = 1/2 (the kissing number case) and 3 < n < 23, n 7^ 8, the Levenshtein 
bounds are better that the Coxeter-Boroczky bounds but weaker than these which 
were obtained by Odlyzko-Sloane [38] . 

In three dimensions, the Levenshtein bound gives T3 < L5(3,1/2) ~ 13.285 
and it can be improved to T 3 < 13.184 which is, of course, not enough. Then 
Anstreicher [3] in 2002 and Musin m in 2003 presented new proofs which were 
based on strengthening the linear programming and using spherical geometry on 
S^. The Musin’s approach will be discussed in more details below. 

In four dimensions, we have T 4 < 7^5(4,1/2) = 26 and this can be improved 
to T 4 < 25.5584 which implies that T 4 is 24 or 25. Then Arestov-Babenko [5] 
proved in 2000 that the last bound is the best possible one can find by pure 
linear programming. Earlier (in 1993), Hsiang |24j claimed a proof that T 4 = 24 
but that proof was not widely recognized as complete. Musin m presented his 
proof of r4 = 24 in 2003 to finally convince the specialists. 

Odlyzko-Sloane [38] use discrete version of the condition (Al) and then apply 
the usual linear programming for s = 1/2 and 3 < n < 24. Their table can be 
seen in m Chapter 1, Table 1.5]. Upper bounds for 25 < n < 32 by linear 
programming were published in m- Now the hrst open case is in dimension 
five, where it is known that 40 < rs < 44 (the story of the upper bounds is: 
T5 < -^5(5,1/2) = 48, Ts < 46.345 from [38]), T5 < 45 from [5] and rs < 44.998 
from [33]). 

Let n and s be fixed, the Levenshtein bound gives A{n, s) < Lm{n, s) and it 
can be improved as seen by Theorem 3. In [8], the first author proposed method 
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for searching improving polynomials f{t) = A^{t)G{t), where A{t) must have 
m + 1 zeros in [—l,s], G{s) = 0 and G{t)/(t — s) is a second or third degree 
polynomial which does not have zeros in [—1, s]. Moreover, one has fi = 0,i =E 
{m, m+1, m+2, m+3} for two or three consecutive coefficients in the Gegenbauer 
expansion of f{t). There restrictions leave several unknown parameters which 
can be found by consideration of the partial derivatives of /(l)//o and numerical 
optimization methods. This approach was realized (see m) by a programme 
SCOD. In fact, SCOD first checks for possible improvements by Theorem 3 and 
then applies the above method. It works well for improving Lm{n, s) for 3 < m < 
16 and wide range of s. 

2.3 Strengthening the linear programming 

The linear programming bounds are based on the following identity 

x,yGC,x^y i=l * j=l VxSC 

where C C is a spherical code, 

k 

i=0 

{vij{x) : j = 1,2,..., r*} is an orthonormal basis of the space Harm(i) of homo¬ 
geneous harmonic polynomials of degree i and = dim Harm(i). In the classical 
case (cf. |16l [25] the sums of the both sides are neglected for polynomials which 
satisfy (Al) and (A2) and this immediately implies Theorem 1. 

Musin 13ZI strengthened the linear programming approach by proposing the 
following extension of Theorem 1 which deals with a careful consideration of the 
left hand side of ([3]). 

Theorem 4. [37] Let f{t) be a real polynomial such that 
(Bl) f{t) < 0 for to <t < s, where —to > s, 

(B2) f{t) is decreasing function in the interval [—l,to]; 

(B3) The coefficients in the Gegenbauer expansion f{t) 
satisfy /o > 0, /fc > 0 /or i = 1,..., m. 

Then 
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where hm, fn = 0 , 1 ,...,/x, is the maximum of f{l) + ei = 

( 1 , 0 ,..., 0 ), over all configurations ofm unit vectors {?/i, 2 / 2 , • • •, Um} the sphe- 
rieal eap (opposite of yi) defined by —1 < {yi,x) < to such that {yi,yj) < s. 

The proof of Theorem 4 follows from ([3|) in a similar way to the proof of 
Theorem 1 - neglect the nonnegative sum in the right hand side and replace the 
sum in the left hand side with its upper bound 

|C| 

j-{yi^yj)<io 

Now observe that the last expression does not exceed ^ 

Now the problems are to find /x, choose to and a polynomial which minimizes 
the maximal value of ho, hi,... ,hfj_. In [37] were found good polynomials f{t) 
by an algorithm which is similar to the algorithm of Odlyzko-Sloane [38]. One 
easily sees that ho = /(I) and hi = /(I) + /(—I). However, the calculation of 
the remaining hm’s usually requires estimations on 

S{n,M) = minmax{s(C') : C C is a spherical code, \C\ = M} 

(cf. P |22l EH [3211371 M\), observe that D{n,M) = ^2(1 - S{n,M))). This 
approach was successfully applied in dimensions three and four. In [37] also noted 
that this generalization does not give better upper bounds on the kissing numbers 
in dimensions 5, 6 , 7 and presumably can lead to improvements in dimensions 9, 
10, 16, 17, 18. 

For n = 3 and s = 1/2 it is provedthat /x = 4, chooses to = —0.5907 and finds 
suitable polynomial of degree 9 (similar to these found in [38l jS] [26]) to show 
that T 3 = 12 . Analogously, for n = 4 and s = 1/2 he has fi = 6, to = —0.608 and 
certain polynomial of degree 9 to obtain T 4 = 24. The calculations of ho, hi and 
/i 2 = maX(^<jr/ 3 {/(l) + / (cos + / (— cos( 7 r /3 — </?))} are easy but computations 
of ho, ■■■ ,ho require numerical methods. 

2.4 Semidefinite programming 

Let C = {xi} C be a spherical code, let I C [—1,1) and let 
SkiC,I):= Y {x,y)^ = \C\Y^tt^ 

{x,y)£l ,x,y£C tGl 

Odlyzko and Sloane [38] used in dimension 17 the constraints 

so{C,Ii)<\C\, so(C,/ 2 ) <2|C|, 
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where Ii = [—1,—\/3/2) and I 2 = [—1, — y^ 2/3), to improve on the LP bound. 
More general, if it is know that the open spherical cap of angular radius (p 
can contain at most m points of C code with S{C) = s, where cosp = t = 
+ (1 — s)/(m + 1), then so{C,I) < m\C\, where I = [—l,t). 

Pfender [39] found the inequality 

S 2 {C,I) <so{C,I)s + \C\{l-s), 


where I = [—1 — y/s), and used it to improve the upper bounds for the kissing 
numbers in dimensions 9, 10, 16, 17, 25 and 26. In fact, the discussion in the 
preceding subsection can be viewed as in the following way: the third author 
[SI ESIEelETI found a few inequalities for some linear combinations of Sk{C,I) 
for 0 < k < 9, s = 1/2 (the kissing numbers’ case) and certain I = [—l,to]) 
to < —1/2. In particular, that gave the proof that T 4 = 24 [37| and a new 
solution of the Thirteen spheres problem |35j . 

This approach can be further generalized by consideration of the three-point 
distance distribution 


Au,v,t = E : {x,y) = u, {x,z) 


V, {y,z) = t}| 


(note that Au^u,i = Au). Here one needs to have 1 -|- 2uvt > 
Bachoc-Vallentin [5] developed this to obtain substantial improvements for the 
kissing numbers in dimensions n = 4,5, 6, 7, 9 and 10. Some numerical difficulties 
prevented Bachoc-Vallentin from furthers calculations but Mittelmann-Vallentin 
|33j were able to overcome this and to report the best known upper bounds so 
far. 


3 Lower bounds on kissing numbers 

3.1 Constructions A and B 

The idea for using error-correcting codes for constructions of good spherical codes 
is natural for at least two reasons - it usually simplifies the description of codes 
and makes easier the calculation of the code parameters. Leech-Sloane |29j make 
systematic description of dense best sphere packings which can be obtained by 
error-correcting codes and give, in particular, the corresponding kissing numbers. 

We describe Constructions A-B following m- Let C be a (n, M, d) binary 
code. Then Construction A uses C to build a sphere packing in R” by taking 
centers of spheres (xi,X 2 , ■ ■ ■, Xn), xi are integers, if and only if the n-tuple 

(xi(mod 2), X 2 (mod 2),... , Xn(mod 2)) 
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belongs to C. 

The largest possible radius of nonoverlapping spheres is ^ min{2, \/d}. The 
touching points on the sphere with center x are 

2‘^Afi{x) if d < 4, 2n + 16^4(x) if d = 4 2n if d > 4, 

where Ai{x) is the numbers of codewords of C at distance i from x. Suitable 
choices of codes for Construction A give good spherical codes for the kissing 
number problem in low dimensions. The record lower bounds for the kissing 
numbers which can be produced by Construction A are shown in Table 1. 

Let in addition all codewords of C have even weight. Construction B takes 
centers {xi^X 2 -, ■ ■ ■ ^Xn), xi are integers, if and only if (xi(mod 2), X 2 (mod 2), 

... ,x,i(mod 2)) G C and 4 divides The touching points on the sphere 

with center x are now 

2 ‘^~^Afi{x) if d < 8, 2n(n — 1) + 128^8(x) if d = 8, 2n(n — 1) if d > 8. 

This, say complication, of Construction A gives good codes for the kissing number 
problem in dimensions below 24. It is remarkable that Construction B produces 
the even part of the Leech lattice in dimension 24. The record achievements of 
Construction B are also indicated in Table 1. 

Having the sphere packings (by Constructions A and B, for example) one can 
take cross-sections to obtain packings in lower dimensions and can build up layers 
for packings in higher dimensions. This approach is systematically used in m 
(see Chapters 5-8). 

3.2 Other constructions 

Dodunekov-Ericson-Zinoviev m proposed a construction which develops the 
ideas from the above subsection by putting some codes at suitable places (sets of 
positions in the original codes; this is called concatenation in the coding theory). 
This construction gives almost all record cardinalities for the kissing numbers in 
dimensions below 24. Ericson-Zinoviev [I8l[lll[20] later proposed more precise 
constructions which give records in dimensions 13 and 14 m- 

4 A table for dimensions n < 32 

The table of Odlyzko and Sloane [T3l|38] covers dimensions n < 24. Lower bounds 
by constructions via error-correcting codes in many higher dimensions can be 
found in [T8l[l3] (see also http://www.research.att.com/~njas/lattices/kiss.html). 
The table below reflects our present (July 2012) knowledge in dimensions n < 24. 


10 


Dimension 

Best known 
lower bound 

Best known 
upper bound 

3 

12 

12 

4 

24 

24 

5 

40 

45 

6 

72 

78 

7 

126 

134 

8 

240 

240 

9 

306 

364 

10 

500 

554 

11 

582 

870 

12 

840 

1357 

13 

1154 

2069 

14 

1606 

3183 

15 

2564 

4866 

16 

4320 

7355 

17 

5346 

11072 

18 

7398 

16572 

19 

10668 

24812 

20 

17400 

36764 

21 

27720 

54584 

22 

49896 

82340 

23 

93150 

124416 

24 

196560 

196560 


The lower bounds in the Table above follow Table 1.5 from m apart from 
dimensions 13 and 14 taken from [20]. The upper bounds are mainly taken from 
[33] (dimensions 5-7, 9-23). 

Note that recently in [12] were found new kissing configurations in 25 through 
31 dimensions, which improve on the records set in 1982 by the laminated lattices. 
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